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5-1 ' Abstract 

In recent years the umbral calculus has emerged from the shadows to provide an elegant correspondence frame- 
work that automatically gives systematic solutions of ubiquitous difference equations — discretized versions of the 
differential cornerstones appearing in most areas of physics and engineering — as maps of well-known continuous 
"T" functions. This correspondence deftly sidesteps the use of more traditional methods to solve these difference equa- 
tions. The umbral framework is discussed and illustrated here, with special attention given to umbral counterparts 
of the Airy, Rummer, and Whittaker equations, and to umbral maps of solitons for the Sine-Gordon, Korteweg-de 
Vries, and Toda systems. 
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^ I Introduction 

p 

t:j- Robust theoretical arguments have established an anticipation of a fundamental minimum measurable 
^ length in Nature, of order L^^^^^^, = yJhG^^/c^ = 1.6162 x 10^'^^ m, the corresponding mass and time being 
Mpj^^^^ = y/hcjG^ = 2.1765 x 10"^ kg and Lpj^^^^/c = 5.3911 x 10"'^'* s. The essence of such arguments is 
the following (in relativistic quantum geometrical units, wherein h, c, and Mp|^_^^^ are all unity). 
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In a system or process characterized by energy E, no lengths smaller than L can be measured, where 
ci ' Lis the larger of either the Schwarzschild horizon radius of the system (~ E) or, for energies smaller than 
the Planck mass, the Compton wavelength of the aggregate process (~ 1/E). Since the minimum of 
max(E, l/£) lies at the Planck mass (E = 1), the smallest measurable distance is widely believed to be of 
order Lpj^^^^j, . Thus, continuum laws in Nature are expected to be deformed, in principle, by modifications 
at that minimum length scale. 

Remarkably, however, if a fundamental spacetime lattice of spacing a = O(Lpj^^^^j) is the structure that 
underlies conventional continuum physics, then it turns out that continuous symmetries, such as Galilei 
or Lorentz invariance, can actually survive unbroken under such a deformation into discreteness, in a 
nonlocal, umbral realization p3llHl- 

Umbral calculus, pioneered by Rota and associates in a combinatorial context \Ty, "41, specifies, in 
principle, how functions of discrete variables in infinite domains provide systematic "shadows" of their 
familiar continuum limit properties. By preserving Leibniz's chain rule, and by providing a discrete 
counterpart of the Heisenberg algebra, observables built from difference operators shadow the Lie alge- 
bras of the standard differential operators of continuum physics. (For a review relevant to physics, see 
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[[loji .) Nevertheless, while the continuous symmetries and Lie algebras of umbrally deformed systems 
might remain identical to their continuum limit, the functions of observables themselves are modified, in 
general, and often drastically so. 

Traditionally, the controlling continuum differential equations of physics are first discretized [J2l[T3}[5|, 
and then those difference equations are solved to yield umbral deformations of the continuum solutions. 
But quite often, routine methods to solve such discrete equations become unwieldy, if not intractible. On 
the other hand, some technical difficulties may be bypassed by directly discretizing the continuum solu- 
tions. That is, through appropriate umbral deformation of the continuum solutions, the corresponding 
discrete difference equations may be automatically solved. However, as illustrated below for the sim- 
plest cases of oscillations and wave propagation, the resulting umbral modifications may present some 
subtleties when it comes to extracting the underlying physics. 

In I'lS] the linearity of the umbral deformation functional was exploited, together with the fact that the 
umbral image of an exponential is also an exponential, albeit with interesting modifications, to discretize 
well-behaved functions occurring in solutions of physical differential equations through their Fourier ex- 
pansion. This discrete shadowing of the Fourier representation functional should thus be of utility in 
inferring wave disturbance propagation in discrete spacetime lattices. We continue to pursue this idea 
here with some explicit examples. We do this in conjunction with the umbral deformation of power se- 
ries, especially those for hypergeometric functions. We compare both Fourier and power series methods 
in some detail to gain further insight into the umbral framework. 

Overall, we utilize essentially all aspects of the elegant umbral calculus to provide systematic solu- 
tions of discretized cornerstone differential equations that are ubiquitous in most areas of physics and 
engineering. We pay particular attention to the umbral counterparts of the Airy, Kummer, and Whittaker 
equations, and their solutions, and to the umbral maps of solitons for the Sine-Gordon, Korteweg-de 
Vries, and Toda systems. 



2 Overview of the umbral correspondence 

For simplicity, consider discrete time, t = 0, a, 2a, • ■ ■ , na, • • ■ . Without loss of generality, broadly 
following the summary review of |To|, consider an umbral deformation defined by the forward difference 
discretization of 9f, 

, , , x(t + a) -x(t) 

Ax t = — , 1 

a 

and whence of the elementary oscillation equation, x{t) = —x(i\ namely, 

l<^x{i) = ^— — ^ ^ ' ^ ' = -x{t) . (2) 

Now consider the solutions of this second-order difference equation. Of course, (dJ can be easily solved 
directly by the textbook Fourier-component Ansatz x{t) a r*, [2], to yield (1 ± ia)'^". However, to 
illustrate instead the powerful systematics of umbral calculus §13} [Toll , we produce and study the solution 
in that framework. 

The umbral framework considers associative chains of operators, generalizing ordinary continuum 
functions by ultimately acting on a translationally-invariant "vacuum" 1, after manipulations to move 
shift operators to the right and have them absorbed by that vacuum, which we indicate by T • 1 = 1. 
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Using the standard Lagrange-Boole shift generator 

7=6"^', so that Tf{t)-l=f{t + a)T-l=f{t + a)l, (3) 

the umbral deformation is then 

T-1 
dt ^ A = , (4) 

a 

t ^ tT-\ (5) 

P I — > [tT-^Y = t{t-a){t-2a)...{t-{n-l)a)T-" = [tYT-'', (6) 

so that [tf = 1, and, for n > 0, [0]" = 0. The [f]" are called "basic polynomials'lll for positive n jT3}[iol |qj, 
and they are eigenf unctions of fT^^A. 

A linear combination of monomials (a power series representation of a function) will thus trans- 
form umbrally to the same linear combination of basic polynomials, with the same series coefficients, 
/(f) I — > f{tT^^). All observables in the discretized world are thus such deformation maps of the 
continuum observables, and evaluation of their direct functional form is in order. Below^, we will be con- 
cluding the correspondence by casually eliminating translation operators at the very end, first through 
operating on the vacuum and then leaving it implicit, so that F [t) = f{tT^^) ■ 1. 

The umbral deformation relies on the respective umbral entities obeying operator combinatorics iden- 
tical to their continuum limit {a — > 0), by virtue of obeying the same Heisenberg commutation relation 



[dt,t] = 1 = [A,tT-^]. (7) 

Thus, e.g., by shift invariance, TAT^^ = A, 

[dt,t"]=nt"-^ ^ [A,[t]"T-"]=n[tY-^T^-", (8) 

so that, ultimately, A[f]" = n[t]"^^. For commutators of associative operators, the umbrally deformed 
Leibniz rule holds HI, 

[AJ{tT-')g{tT-')] = [A,f{tT-'MtT-')+f{tT-')[\g{tT-')] , (9) 

ultimately to be dotted onto 1 . Formally, the umbral deformation reflects (unitary) equivalences of the 
unitary irreducible representation of the Heisenberg-Weyl group, provided for by the Stone-von Neu- 
mann theorem. Here, these equivalences reflect the alternate consistent realizations of all continuum 
physics structures through systematic maps such as the one we have chosen. It is worth stressing that 
the representations of this algebraic relation on the real or complex number fields can only be infinite 
dimensional, that is, the lattices covered must be infinite. 

Now note that, in this case the basic polynomials [t]" are just scaled falling factorials, for n > 0, i.e. 
generalized Pochhammer symbols, which may be expressed in various w^ays: 



[t]" = (tT-^Y ■l = t{t-a)---{t-{n-l)a) 



{t/a-n)\ 

Thus [-t]" = {-)"[t + a{n-l)Y. Furthermore, [an]" = a"nl ; [tY[t - am]"-"" = [t]" for < m < n ; and 
for integers < m < n, [am]" = 0. Thus, A'"[f]" = [anYlt]"-"" /a"'. 



^We stress that the notation [t]" is shorthand for the product t{t ~ a)...{t — {n — l)fl). It is not just the wth power of [t] = t. 
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Negative umbral powers, by contrast, are the inverse of rising factorials, instead: 

- ■t-\n / -, \ n 



t 



= irlV-i=, .. ' . . = »-» ('/"'i 



a 



tj {t + a){t + 2a)--- {t + na) {t/a + n)\ 

r(i + f. + i) ^ ' r(-i) 

These correspond to the negative eigenvalues of tT^^A. 

The standard umbral exponential is then natural to define as |T3| [6|H 

E(Af,Aa) ^ .^W ^ .^^^" . 1 = £ ^[tr = ti^^rCi") = (1 + Aa)'/^ (12) 

11=0 "• n=o V " / 

the compound interest formula, with the proper continuum limit (a — > 0). N.B. There is always a o at 
A = -1/fl. 

Evidently, since A ■ 1 = 0, 

Ae^W = A /W, (13) 

and, as already indicated, one could have solved this equation directljo to produce the above E{At, Aa). 
Serviceably, the umbral exponential E happens to be an ordinary exponential, 

and it actually serves as the generating function of the umbral basic polynomials, 

a" 



^(l + ^«)^'" 



9A 



= W- (15) 

A=0 



Conversely, then, this construction may be reversed, by first solving directly for the umbral eigenfunction 
of A, and effectively defining the umbral basic polynomials through the above parametric derivatives, in 
situations w^here these might be more involved, as in the next section. 

As a consequence of linearity, the umbral deformation of a power series representation of a function 
is given formally by 



/(f) ^ f(f)^/(fT-i)-l = /(A) (l + Afl)^/« 



(16) 

A=0 



This may not alw^ays be easy to evaluate, but, in fact, the same argument may be applied to linear 
combinations of exponentials, and hence the entire Fourier representation functional, to obtain 



m-j_Jrmj_^-e-^il^rc.ar^^ (l + «-) /(. 



(17) 



The rightmost equation follows by converting ico into 9t- derivatives and integrating by parts away from 
the resulting delta function. Naturally, it identifies with Eqn ( fi6l ) by the (Fourier) identity f{^x)g{^) \x=q = 



^Again we stress that e'^M is a short-hand notation, and not just the usual exponential of A [t] = \t. 

3N.B. There is an infinity of "non-umhral" extensions of the E {M,\a) solution [9]: Multiplying the umbral exponential by 
an arbitrary periodic function g{t + a) = g{t) will pass undetected through A, and thus will also yield an eigenfunction of A. 
Often, such extra solutions have either a vanishing continuum limit, or else an ill-defined one. 
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g{dx)f{x) \x=o- It is up to individual ingenuity to utilize the form best suited to the particular application 
at hand. 

It is also straightforward to check that this umbral transform functional yields 

dtf ^ Af , (18) 

and to evaluate the umbral transform of the Dirac delta function, which amounts to a cardinal sine or 
sampling function, 

smj^jl + t/a)) 
'^'^ ^ (nia + t)) ' (^9) 

or to evaluate umbral transforms of rational functions, such as 

/ = T^^ ^ F = e^l'a"''Y(tla + 1, l/a) , (20) 

to obtain an incomplete Gamma function (A&S 6.5.3), ^^id so on. Note how the last of these is distinctly, 
if subtly, different from the umbral transform of negative powers, as given in (jTT)) . 

In practical applications, evaluation of umbral transforms of arbitrary functions of observables may 
be more direct, at the level of solutions, through this deforming functional, Eqn ( [17] ). For example, one 
may evaluate in this way the umbral correspondents of trigonometric functions, 

Sin[f] ^ ^. , Cos[i] ^ Y , (21) 

so that 

ASin[i] = Cos[f] , ACos[f] = -Sin[i] . (22) 

As an illustration, consider phase-space rotations of the oscillator. The umbral deformation of phase- 
space rotations, 

x = ^, p = -x ^ AX{t)=P{t), AP{t) = -X{t), (23) 

readily yields, by directly deforming continuum solutions, the oscillatory solutions, 

X{t) = X(0) Cos[f] + P(0) Sin[f], P{t) = P(0) Cos[f] - X(0) Sin[f] . (24) 

In view of ( [14] ), and also 

(1 + ia) = a/TT^ e<arctan(«) ^ (25) 

the umbral sines and cosines in ( [24I are seen to amount to discrete phase-space spirals, 

X{t) = (1 + a^)^(x{0) cos{cvt) + P(0) sin(a;f)) , P{t) = (1 + a^)^ (p{0) cos{cot) - X(0) sin(a;f)) , 

(26) 
w^ith a frequency decreased from the continuum value (i.e. 1) to 

CO = arctan(fl)/fl < 1 . (27) 

So the frequency has become, effectively, the inverse of the cardinal tangent function|f| Note that the 
umbrally conserved quantity is, 

2S = X(0)2 + P(0)2 = (1 + a^)^(x{tf + P{tf) , (28) 



such that AS = 0, with the proper energy as the continuum limit. 



4That is, for = arctan(fl), the spacing of the zeros, period, etc, are scaled up by a factor of tanc(0) = ^']k" > 1 • For 
complete periodicity on the time lattice, one further needs return to the origin in an integral number of N steps, thus a solution 
oi N = 2nn/ arctana. Example: For a = 1, the solutions' radius spirals out as 2'^^, while cv = n/A, and the period is t = 8. 
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3 Reduction from second-order differences to single term recursions 

In this section and the following, to conform to prevalent conventions, the umbral variable will be de- 
noted by X, instead of t. In this case there is a natural way to think of the umbral correspondence that 
draws on familiar quantum mechanics language [13] : The discrete difference equations begin as operator 
statements, for operator xs and Ts, but are then reduced to equations involving classical-valued func- 
tions just by taking the matrix element (^| • ■ ■ \vac) where \vac) is translationally invariant. The overall 
x-independent non-zero constant {x\vac) is then ignored. 



To be specific, consider Whittaker's equation (A&S 13.1.31) for ^ = 1/2, 



9x + ^-J)yW = o. (29) 

This umbrally maps to the operator statement 

A' + T^-l)y(xT-^)=0. (30) 

Considering either y{xT~^) ■ 1 = Y{x), or else {x\ y{xT^^) \vac) = y{x) {x\vac), this operator statement 
reduces to a classical difference equation, 

2 2 

Y{x + 2a) - 2Y{x + a) + Y{x) + -^^ Y{x + a) - — Y{x) = . (31) 

Before using umbral mapping to convert continuous solutions of (J29]) into discrete solutions [TT.TT] of 
(J31II, here we note a simplification of the latter equation upon choosing a = 2, which amounts to setting 
the scale of x. With this choice ( [31] ) collapses to a mere one-term recursion. Shifting x — )■ x — 2 this is 

y(x + 2)=2(^-^)y(x). (32) 

Despite being a first-order difference equation, however, the solutions of this equation still involve two 
independent "constants of summation" even for x restricted to only integer values, because the choice 
fl = 2 has decoupled adjacent vertical strips of unit width on the complex x plane. To be explicit, for 
integer x > 0, forward iteration gives [2I 



^L_^'|y(l) and Y {2k + 2) = 2'^ (j\ ^ 



y(2/: + l) =2Mn^— ^ ) y(l) and y(2^ + 2) = 2M fj^^ ) y (2) , for integer fc > , (33) 



with y (1) and Y (2) the two independent constants that determine values of Y for all larger odd and even 
integer points, respectively. 

Or, if generic x are contemplated, the equation ([32]) has elementary solutions, for arbitrary complex 
constants Ci and C2, given by 

2x/2-p/'x_^^ (_2Y^^ 



r(j) \ ■ n' ' - \2 
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In the second expression, we have used T (z) F (1 — z) = tz/ sin nz. Note the C2 part of this elementary 
solution differs from the Ci part just through multiplication by a particular complex function with period 
2. This is typical of solutions to difference equations since any such periodic factors are transparent to A, 
as mentioned in an earlier footnote M. 

As expected, even for generic x the constants Ci and C2 may be determined given Y (x) at two judi- 
ciously chosen points, not necessarily differing by an integer. For example, if < k < 1, 



Ci = ^^Y(2 + 2K), C2 



^^^Ci-^F(k) y(2). 
sinTTK 2 ^ ^ ^ ^ 



(36) 



Moreover, poles and zeros of the solution are manifest either from the F functions in ( [34] ), or else from 
continued product representations such as ([33]). For the latter, either forward or backward iterations of 
the first-order difference equation ([32]) may be used. Schematically, 



or alternatively. 



y(x) 



Y(x) 



(2j: - 4 - 4k) (2x - 8 - 4k) {2x - 12 - 4k) 
(x-2)(x-4)(x-6)--- 

x{x + 2){x + 4:) ■■ ■ 
' (2x - 4k)(2;c + 4(1 - k)){2x + 4(2 - k)) ■ 



(37) 



(38) 



Although both terms in ( [34] ) have zeroes, the Ci term also has poles while the C2 term has none — it 
is an entire function of x — and it is complex for any nonzero choice of C2. Of course, since the equation 
( [32] ) is linear, real and imaginary parts may be taken as separate real solutions. All this is evident in the 
following plots for various selected integer k. 




r(^-) 



for K = 1, 2, and 3 in red, blue, and green. 
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2'^^ cos tt{^x) 

r(ix)r(i-ix+)c) 



for K = 1, 2, and 3 in red, blue, and green. 



4-- 



-2 -1 



-3 




rfi \rft 1 - \ ^°^ K^ = 1a 2, and 3 m red, blue, and g 



reen. 
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Collapse to a mere one-term recursion also occurs for an inverse-square potential. 



K 



^' + ^ - ^j l/(^) = ■ (39) 

For /^fl^ = 1, which amounts to setting the scale of the energy of the solution, the umbral version of this 
equation reduces to 

That is to say. 

Elementary solutions for generic x, for arbitrary complex constants Ci and C2, are given by 



Again, the C2 part of this elementary solution differs from the Ci part just through multiplication by a 
particular complex function with period a. And again, poles and zeros of these and other solutions are 
manifest either from those of the T functions, or else from a continued product form, e.g. 

, . _ {x^ + xa + Ka^) {{x + a)^ + {x + a)a + Ka^)- ■■ 

^''' ~ {2x{x + a)) {2{x + a)ix + 2a))--- " ^44) 

It is not surprising that ([29]) and ( [39] ) share the privilege to become only first-order difference equations 
for specific choices of a, as in ([32]) and ([41]), because they are both special cases of Whittaker's differential 



equation, as discussed in the next section. No other linear second-order ODEs lead to umbral equations 
with this property. 



4 Discretization through hypergeometric recursion 

In this section we discuss several examples using umbral transform methods to convert solutions of 
continuum differential equations directly into solutions of the corresponding discretized equations. We 
use both Fourier and power series umbral transforms. 

As an explicit illustration of the umbral transform functional fi7| ), inserting the Fourier representation 



of the Airy function (A&S 10.4.32) yields 



AiryAi (x) 1 — > UmAiryAi {x, a) =Re(— e^'*^' (1 + jka) » dfc j . (45) 



This integral is expressed in terms of hypergeometric functions and evaluated numerically in Appendix A 
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Likewise, gaussians also map to hypergeometric functions, as may be obtained by formal series ma- 
nipulations: 



G{x,a) = Y^ 



)"[xF" ~ 1 



n=0 



nl 



Y^±.^-iYa^n r(f + i) 



n=Q 



n\ 



r 



2n + l) 



gr(n-if)r(n + i-if: 



-4fl 



2\n 



n=0 



r(-^f) r( 



1 _ 1 x^ 

2 la, 



2^0 



2 fl 2 V a 



where the reflection and duplication formulas were used to write 



r(^ + i) _4"r(n-if)r(n + i-if] 



r(f-2n + l^ 



r(-i§)ra-if) 



(46) 

(47) 
(48) 

(49) 



la, 



.2 2/ 



While the series ( [47] ) actually has zero radius of convergence, it is Borel summable, and the resulting 
regularized hypergeometric function is well-defined. See Appendix B for some related numerics. 



For another example drawn from the familiar repertoire of continuum physics, consider the confluent 

(50) 



hypergeometric equation of Kummer (A&S 13.1.1) 



xy" + {^-x) y' -ay = Q, 



whose regular solution at x = 0, jjexpressed in various dialects is 



y = iFi (a; /5; x) = M (a, /3, x) = KummerM (a, /5, x) , 
with series and integral representations 

.fr^ r (a) r(/3 + n)n! 



(51) 



(52) 



n=0 



r(/3) 



r(a)r(/5-a) 



^ + /3^+2/3(/3 + l) 



1 ^(^ + 1) ^2 , l ^(a + l)(« + 2) 3 / 4 

2/3(6 + 1) ^6/3(/3 + l)(/3 + 2) ^ "^ ^ 



The second, independent solution of ( [go] !, with branch point at x = 0, is given by Tricomi's confluent 
hypergeometric function (A&S 13.1. 3). sometimes known as HypergeometricU^: 



U{oi,^,x) 



n 



M{a,^,x) 



sin TTjB vr(l + a-|6)r(^) 



i_^ M(l + ^-/3,2-/3,x) 

r(«)r(2-/3) 



(53) 



Invoking the umbral calculus for x, either of these confluent hypergeometric functions can be mapped 
onto their umbral counterparts using 



iFi (a;j6;x) 



2F1 (a, — ;,6;-fl) , 



(54) 



where 2F1 is the well-known Gauss hypergeometric function (A &S 15.1.1) , This map from ifi to 2^1 
follows from the basic monomial umbral map. 



xT- 



T^ + l 






(55) 
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and from the series ( [g^] ). When combined, these give the well-kno^vn series representation of 2^1 • 

Next, reconsider the one-dimensional Coulomb problem defined by Whittaker's equation for general 



ji (A&S 13.1.31) 



//, ( 1 , ^ , il-l^'' 



y"+^-4 + - + -^^^Jy = o. (56) 

Since k and ]i are both arbitrary, this also encompasses the inverse-square potential, ( [39] ). lExact solutions! 
of this differential equation are 

y (x) = Ci whittakerM (k, ji, x) + C2 whittakerW (k, /^, x) , (57) 



whittakerM (k, pi, x) = a:f'+i/2g-x/2 ^^^ / ^ _ ^ + _; 2/^ + 1; x j , 



/ FrlrSiTifi(?^-'^ + 2;2/^ + i;^) \ 



whittakerW (K,/i,x) = x 



}i+l/2 -x/2 






(58) 



(59) 



V +F(^ ^"'' ^^1 (-?^ - '^ + ^; -2f^ + 1;-) y 



Umbral versions of these solutions are complicated by the exponential and overall power factors in the 
classical relations between the ifi's and the Whittaker functions, but this complication is manageable. 
(In part this is because in the umbral calculus there are no ordering ambiguities |T7|.) 

To obtain the umbral version of the Whittaker functions, we begin by evaluating 

CO 00 / -1 \ m Hiill) r -im+n 
g-ixT-i r (...a.^^-l\ , - V- V- I 1\ r(«) m 



cooo/-|\m \ / 
m=On=0 V / Y(R) 



m\n\ 

a\-a ^ ( X ^ la . ,^ ^ 



where we have performed the sum over m first, to obtain 



1\"' 1 fl'" 1 / a^ ^-"" 



£oV 2J r(f-n-m + l)m! r(f-n + l)(^ 2) " ^^'^ 

The sum over n then gives the Gauss hypergeometric function in ( [6ol) . 

Next, to deal with the umbral deformations of the Whittaker functions, we need to use the continuation 
of (tToJ and (jTTj) to an arbitrary power of xT^^, namely. 

This continuation leads to the following: 



V a 



^-i(.-7«)T-^ iFi (a;/3; (x - 7fl) T-') ■ 1 . (63) 
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Thus we obtain the umbral map 

A-*',F,(«;^;.)^^l(i±iL_«.(l-|)'-%f, {^,7-'-:P.^) . (64) 

Finally then, specializing to the relevant oc, ji, and 7, we find the umbral Whittaker functions. In particu- 
lar, 

whittakerM(K,;^,x)^— ^^^^-fl^+i/2(^l--j 2F1 (^ + 2 " '''^ + 2 " «'^ + ' ^^j ' 

(65) 

This result for general a exhibits what is special about the choice fl = 2, as exploited in the previous 
section. To realize that choice from ^^ requires taking a limit a /^ 2, hence it requires the asymptotic 



behavior of the Gauss hypergeometric function (A&S 15.3.7) 

2^-00 r (p) r (7 — a) r (a) r (7 — j6) 



Now with sufficient care, a = 2 solutions can be coaxed from the umbral version of whittakerM in ( [65) , 
and /or the corresponding umbral counterpart of whittaker W, upon taking \m\i,/2 arid making use of 
(|66l >. Moreover, in principle the umbral correspondents of both Whittaker functions could be used to 
obtain from this limit a solution with two arbitrary constants. 

On the other hand, for a = 2, the umbral equation corresponding to ([56]) again reduces to a one-term 
recursion, namely, 

Y(,c + 2) - 2(x + 2)(x-2k) 

For generic x, solutions for arbitrary complex constants Ci and C2 are then given by 



r(f + 5 + j')r(f + J-F) V t' ^2^ ^2 



w^hich agrees with ( [34] ) when }i = 1/2, of course. As in that previous special case, the C2 part of (|68l > 
differs from the Ci part just through multiplication by a particular complex function with period 2 IJ9J. 

We graph some examples to show the differences between the Whittaker functions and their umbral 
counterparts, for a = 1. 
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whittakerM (k, 1/2, x) for k = 1, 2, and 3 in red, blue, and green. 




Umbral whittakerM (k, 1/2, x) for a = 1, and for k = 1, 2, and 3 in red, blue, and green. 
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The examples above are specific illustrations of combinatorics that may be summarized in a few umbral 
hypergeometric mapping lemmata, the simplest being 



Lemma i: 



pFq (ai,...,ap;,6i,...,,6q;x) 1 — > p+iFc, {oci,...,oip,-x/a;^i,...,Pq;-a) , (70) 

where the series representation of the generalized hypergeometric function pfg i^ 

,F,[a„ ,^,,fi„ 'P^''^-T{a,)...Y{a,)tonfii+n)---r{fi, + n)nr ^^D 

A proof of ( [70] ) follows from formal manipulations of these series. 

The umbral version of a more general class of functions is obtained by replacing x — )■ xT^^ in functions 
of x^ for some fixed positive integer k. Thus, again for hypergeometric functions, we have 



Lemma 2: 



pFci[ocir-- ,cip]^ir-- ,^q;x^j 1 — > (72) 

And again, a proof follows from formal series expansions. 

Multiplication by exponentials produces only minor modifications of these general results, as was 
discussed above in the context of Whittaker functions, namely. 



Lemma 3: 



s^%F<, (*i,--- ,ap;i6i,--- ,i6,;x'^j ^ (73) 

1 / x\ 1 /^ x\ 1 /, . x\ „ „ / —ak ^ 



In addition, multiplication by an overall power of x gives 



Lemma 4: 



x^e^"" pF^l^di,--- ,(x.p;^ir-- ,Pq;x^) i — > (74) 

Y{^ + l)a^{l + a\y-^ ( 7-f 1 + 7-f ^_zl±l^. a R ■ ( ''^ \'\ 

r(i-7 + l) P+^^? (^«i'---'V ^ ' ^ '■■•' ^ 'P^'"' 'P'^'\l + ax) ) ■ 

5Recall results from using the ratio test to determine the radius of convergence for the pF,j (aj, • • • , oc-p; /3i, ■ ■ ■ , /5,j; x) series: 
li p < q + 1 then the ratio of coefficients tends to zero. This implies that the series converges for any finite value of x. 
li p = q + 1 then the ratio of coefficients tends to one, hence the series converges for jxj < 1 and diverges for |x| > 1. 
It p > q + 1 then the ratio of coefficients grows without bound. The series is then divergent or asymptotic, and is a symbolic 
shorthand for the solution to a differential equation. 
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5 Wave propagation 

Given the umbral features of discrete time and space equations discussed above, separately, it is natural 
to combine the two. 

For example, the umbral version of simple plane waves in 1+1 spacetime would obey an equation of 
the type HI [9), 

(a2-A?)F = 0, (75) 

on a time-lattice with spacing a and a space-lattice with spacing b, not necessarily such that b = a in all 
spacetime regions. For generic frequency, wavenumber and velocity, the basic solutions are 

For right-moving waves, say, these have phase velocity 

, ,, CO flarcsinffc) , , 

v{cv,k) = - r^ . (77) 

k e'arcsm(fl) 

Thus, the effective index of refraction in the discrete medium is (foarcsin(fl))/(flarcsin(b)), i.e. modified 
from 1. Small inhomogeneities of a and b in the fabric of spacetime over large regions could therefore 
yield interesting effects. 

Technically, a more challenging application of umbral methods involves nonlinear, solitonic phenom- 
ena I118J , such as the one-soliton solution of the continuum ^Sine-Gordon equation^ 



(di - df)f {x, t) = sin(/ (x, 0) , /SG (x, t) = 4arctan (meV^n . (78) 



The corresponding umbral deformation of the PDF itself would now also involve a deformed potential 
sm{f{xT^^,tTfr^)) ■ 1. But rather than tackling this difficult nonlinear difference equation, one may 
instead use the umbral transform ( [17] ) to infer that /sg {x, t) maps to 

FsGia,b;x,t)= T^^^I^ 4arctan fme^) e'^'^-'^'^l + icoaY^' (1 - ikb)'^^ . (79) 

J-00 (27r) V / 



The continuum Korteweg-de Vries soliton is likewise mapped: 



fKdv{x,t) = -sech { — {x-vt)) 1 — > (80) 

FKdvKfc;.,0=£^^^^ech^(f U-.T)) e-^-^il + icoaf^il-rkbr^' . 

Closed-form evaluations of these Fourier integrals are not available, but the physical effects of the dis- 
cretization could be investigated numerically, and compared to the Lax pair integrability machinery of 
[[lojl , or to the results on a variety of discrete KdVs in [ [14] , or to other studies [7]. 

However, a more accessible example of umbral effects on solitons may be found in the original 
[Toda lattice model I.i6:l . For this model the spatial variable is already discrete, usually with spacing 
& = 1 so X = n is an integer, while the time t is continuous. The equations of motion in that case are 

dq{n,t) dp{n,t 



dt 



p („, t) , ^^^ = - (e-('?("+i'')-''("'^)) - e-(iKO-?("-i.O)) , (81) 
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for integer n. Though x = n is discrete, nevertheless there are exact multi-soliton solutions valid for all 
continuous t, as is well-known. 

Specific one-soliton Toda solutions are given for constant oc, fi, 7, and qo by 

/ 1 + aexp ( — /3n + 7O \ ,^ ^ 

q (n, t) = qo + log / ,/ ,/ ^ , , 82 

( e-n^-rt g-(n+l)/3+7t \ 



■provided that 



7 = ±2sinh(|) . (84) 



So the soliton's velocity is just y = ± 4 sinh f ^ j . 

While obtained only for discrete x = n, for plotting purposes q [n, t) may be interpolated for any x 
(see graph below). To carry out the complete umbral deformation of this system, it is then only necessary 
to discretize t in the equations of motion ((81)) . Consider what effects this approach to discrete time has 
on the specified one-soliton solutions. 

To that end, expand the exact solutions in JSil) as series, 

q (n, = ^0 + E ^ 1—^ (^"'^ - 1) exp {ikt) . (85) 

Upon umbralizing t, the one-soliton solutions then map as 

q (n, ^ Q {n, = ^0 + E ^ IT^ ( ^"'^ " 1 ) (1 + ikay/" , (86) 

k=i '^ ^ ^ 

and these are guaranteed to give solutions to the umbral operator equations of motion, 

Aq (n, f r-i) = 1 (T - 1) ^ (n, fT-^) = p (n, tT'^^ , (87) 

Ap (n, fT-i) = i (T - 1) p (n, tT'^) = - (^e-i<"+'''^'')-<''''^'')) - e-i^"'^^-^)^"-''^^-'))^ , (88) 
upon projecting onto a translationally invariant "vacuum" (i.e. Q [n, t) = q [n, tT^^^ ■ 1). 
Now, for integer time steps, t/a = m, consider the series at hand: 

S(m,c,z) = E|(e-^^-l)(l + cfc)'"=ln(^-l^)+Ec^('j)R(;,z), (89) 

where c = ya, z = —oce~^^, and where for j > 0, 

CO 

R {h z) = E (^"'^ - l) z'^^"' = ^ {''^^' 1 - ;. O) - O (z, 1 - /, 0) . (90) 

jt=0 

Fortunately, for positive integer t/a, we only need the ILerch transcendent functionl 

oo k 
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for those cases where the sums are expressible as elementary functions. For example. 



)t=o ^~^ k^o (l-z) )t=o (1-z) fc=o (1 



-r 



(92) 



The In ( ■ ■ • ) term on the RHS of ([89]) then reproduces the specified classical one-soliton solutions at f = 0, 
while the remaining terms give umbral modifications for t j^ 0. 



Altogether then, we have 



Q{n,t = ma) = g (n, 0) + £ {-fay h) ( O (-ae-/'(«+i), 1 _ y, 0) - O (-ae"^", 1 - ;, 



(93) 



These umbral results are compared to some time-continuum soliton profiles for t/a = 0, 1, 2, 3, and 4 in 
the following Figure (with q^ = 0, a = 1 = ^, and 7 = 2sinh (1/2) = 1.042 ). 



q and Q 




4 5 



Toda soliton profiles q interpolated for all x G [—5, 5] at integer time slices superimposed with their time 

umbral maps Q (thicker curves) for a = 1. 



Thus, the umbral-mapped solutions no longer evolve just by translating the profile shape. Rather, 
they develop oscillations about the classical fronts that dramatically increase with time, that evince not 
only dispersion but also generation of harmonics, and that, strictly speaking, disqualify use of the term 
Isolitonl for their description. Be that as it may, this model is referred to in some studies as integrable ^. 



Curtright &l Zachos Umbral Vade Mecum l8 

These umbral effects on wave propagation evoke scattering and diffraction by crystals. But here the 
"crystal" is spacetime itself. It is tempting to speculate based on this analogy. In particular, were a well- 
formed wave packet to pass through a localized region of crystalline spacetime, with sufficiently large 
lattice spacings, the packet could undergo dramatic deformations in shape, wavelength, and frequency — 
far greater than and very different from what would be expected just from the dispersion of a free packet 
propagating through continuous space and time. 



6 Concluding remarks 

We have emphasized how the umbral calculus has visibly emerged to provide an elegant correspondence 
framework that automatically gives solutions of ubiquitous difference equations as maps of well-known 
continuous functions. This correspondence systematically sidesteps the use of more traditional methods 
to solve these difference equations. 

We have used the umbral calculus framework to provide solutions to discretized versions of several 
differential equations that are widespread building-blocks in many if not all areas of physics and en- 
gineering, thereby avoiding the rather unwieldy frontal assaults often engaged to solve such discrete 
equations directly. 

We have paid special attention to the Airy, Kummer, and Whittaker equations, and illustrated several 
basic principles that transform their continuum solutions to umbral versions through the use of hyperge- 
ometric function maps. The continuum limits thereof are then manifest. 

Finally, we have applied the solution-mapping technique to single solitons of the Sine-Gordon, Korteweg- 
de Vries, and Toda systems, and we have noted how their umbral counterparts — particular solutions of 
corresponding discretized equations — evince dispersion and other non-solitonic behavior, in general. 
Such corrections to the continuum result may end up revealing discrete spacetime structure in astrophys- 
ical wave propagation settings. 

We expect to witness several applications of the framework discussed and illustrated here. 
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Appendix A: Umbral Airy functions 

Formally, these can be obtained by expressing the Airy functions in terms of hypergeometric functions 
and then umbral mapping the series. The continuum problem is given by 



where 



y" - xy = , y (x) = d AiryAi (x) + Ci AiryBi (x) , (94) 

AiryAi (X) = 3,,3r'2/3) '^' (' I V") " WWi^ ''' (' t' l'') ' ^^5) 

AiryBi (X) = 3,,,r'(2/3) '^' (' I V") + T^t) °^^ (' ^ ^"') ' ^''^ 

The y i — )■ Y umbral images of these, solving the umbral discrete difference equation [91 [3| 

y(x + la) - 2Y(x + fl) + y(x) - a^xY[x + a) = , (97) 

are then given by ( J72I 1 for /c = 3. In particular, 

UmAiryAi (x,«) = ^,73^^ 3^1 (-^^, ^ (^ " 3 4 (^ " D ' ^' "^^^ 

3fi -:;-,;; 1--,:; 2-- ;-;-3fl3 . (98) 



3i/3r(l/3) -" ' V 3fl'3 V a)'?>\ aJ'?> 

Since the number of "numerator parameters" in the hypergeometric function ^Fi exceeds the number of 
"denominator parameters" by 2, the series expansion is at best asymptotic. However, the series is Borel 
summable. In this respect, the situation is the same as for the umbral gaussian (see Appendix B). 

Alternatively, as previously mentioned in the text, using the familiar integral representation of AiryAi (x), 
the umbral map devolves to that of an exponential. That is to say. 



AiryAi (xT^M = ^ exp f -is^ + isxT^'^ 



ds (99) 



1 /■+~ 1-3 I 

UmAiryAi (x,fl) = -— / ga'^ (l + isa)" ds . (100) 

277 J — CO 

Just as AiryAi (x) is a real function for real x, UmAiryAi (x, a) is a real function for real x and a, 

( 1 /'+~ 1-3 1- \ 

UmAiryAi (x, a) = Re I — / es'^ (1 + isa) " ds\ . (101) 

After some hand-crafting, the final result may be expressed in terms of just three 2^2 generalized 
hypergeometric functions. To wit. 

Re I- /^~ e¥'^ (1 + isw)-^ d^ (102) 

^, . ( 8a;2sin(7rz/3) Ci(w,z)Hi (w,z) -12z(;(l + 2cos(27rz/3)) C2(z«,z)H2(zi;,z) \ 
^^"^'"''^H +3C3(r.,z)H3(z<;,z) )' 
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where the hypergeometric functions 2^2 («/ b; c, d; z) appear in the expression as 

/1\/1 1 \ fl 1 1121\ 

Hi (z.,z) = r (^-zj r ^- + -zj 2F2 ^-z, - + -z; -, -; ^J , (103) 

H3Kz)=r0 + iz)r(l + iz)2F2g + ^z,l + lz;i4;3i3), (X05) 

and where the coefficients in ( I102I) are 

1 e-(^ln3+llnz<;2)z 

^° ^"'' ^^ " 96^ (sin Ittz) F (Iz) (sin l/r (1 + z)) F (1 + Iz) (cos i/r (2z + 1)) F (| + Iz) ' ^^°^^ 



Ci (w,z) = 2v^cos I -TTZ — -zsignum {w) n j + v^cos i -nz+ -zsignum (zv) n 

+ 36 sm I -TTZ + -zsignum (if) TT I , (107) 

(111 \ /111 \ 

-TT + -TTZ + -TTZ signum (zf) I + v^sin I -n + -ttz + -Trzsignum (iv) | , (108) 



mum zv 



C3 (iti,z) = 3v3cos I -TTZ— -TTZ signum (w) j — 6-^3 cos I -7rz+ -Trzsig 

+ 6-^3 cos I -TTZ + -TTZ signum (w) j + 3 x 33 sin ( -ttz — -Trzsignum {w) j 

— 3v^cos i-KZ — -Trzsignum (if) j + 2 x 3' sin | -ttz + -zti signum (if) j 

+ 3 3 sm I - TTZ H — TTZ signum [w) j +2x33 sm ( - ttz H — ttz signum (if) j . (109) 

While the coefficient functions C0-3 are not pretty, they are comprised of elementary functions, and they 
are nonsingular functions of z. On the other hand, the hypergeometric functions do have singularities 
and discontinuities for negative z. However, the net result for UmAiryAi is reasonably well-behaved. 

We plot UmAiryAi {x,a) for a = 0, ±\, ±\, and ±1. 
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Umbral Airy 




UmAiryAi {x,a) for a = ±1, ±1/2, and ±1/4 (red, purple, & blue dashed/solid curves, resp.) compared 

to AiryAi (x) = UmAiryAi (x, 0) (black curve). 
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Appendix B: Umbral gaussians 

As discussed in the text, straightforward discretization of the series yields the umbral gaussian map: 

~ (-)"M^" 



G(x,fl) = ^ 

n=0 



n\ 



oo ( _'\n 

E ^:t- x{x-a)---{x- {In - 1) a) 



n=0 



n\ 



2^0 



If 1 /"l _ f 

la'lX a 



-Aa' 



(110) 

(111) 



(NB G {x,a) 7^ G (— x, a).) Now, it is clear that term by term the series (|iio|) reduces back to the 
continuum gaussian as fl — ?► 0. Nonetheless, since the series is asymptotic and not convergent for \a\ > 0, 
it is interesting to see how this limit is obtained from other representations of the hypergeometric function 
in (|iii|) , in particular from using readily available numerical routines to evaluate 2^0 for specific small 
values of a. Some examples are shown here. 




G [x, 1/2") versus x S [—3,2], for n = 1, 2, and 3, in red, blue, and green, respectively, compared to 

G {x,0) = exp (— ^^), in black. 



Mathematica® code to produce these graphs is available online, for those interested. It is amusing 
that Mathematica manipulates the Borel regularized sum to render the 2^0 in question in terms of Tri- 
comi's confluent hypergeometric function U, as discussed above in the context of Kummer's equation, cf. 
( [53] ). Thus G can also be expressed in terms of iFi s. The relevant identities are: 

-»)'2-4?j = ^^(^«) [ r(i(i-f)) J?HS )■ 

(112) 



G (x,fl) 



(2fl)«"^LZ 
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